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Abstract. We consider the Tsallis distribution as a source of the apparent 
slope of one-particle spectra in heavy-ion collisions and investigate the equation 
of state of this special kind of quark matter in the framework of non-extensive 
thermodynamics. We relate the energy per particle to the power-law tail of 
spectra at a given temperature. 
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1. Introduction 

1.1. Why temperature 

There is a longstanding discussion in heavy-ion physics whether the concept of 
temperature and the assumption of thermal equilibrium can be applied for the 
description of the fireball formed before hadronisation in high-energy reactions. We 
address here this question from the viewpoint of very rough bulk properties relating 
the inverse slope of one-particle transverse momentum spectra and the energy per 
particle. Hadron mixture fits, the so called "thermal model" assumes a constant 
temperature all over the fireball and obtains the energy per particle using integrals 
with thermal distributions: basically the Gibbs distribution, which leads to Fermi 
and Bose distributions in the grandcanonical treatment. Although it has been 
shown that due to the finiteness of the available phase space volume a canonical 
suppression may occur the existence of temperature was not questioned in this 
framework. 



Study of the one-particle spectra gives a reason to believe in the temperature as 
a universal parameter of the Gibbs distribution, exp(—E/T). From AGS through 
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SPS to RHIC experiments the transverse momentum spectra of several particles, 
first of all pions, kaons and protons, show an exponential fall. This occurs in its 
purest form in the pt — 1 — 4 GeV range. At lower momenta deviations can be 
seen, traditionally attributed to the presence of a global transverse flow. We do not 
intend to discuss this part of the spectra here. At higher momenta in the RHIC 
experiment, however, deviations occur too. Instead of the exponential curves a 
power-law tail fits the experimental hadron data in the p t = 4 — 12 GeV range. 

1.2. Why no temperature 

The power law tail in the pion spectrum contradicts to the assumption of Gibbs 
distribution. Although one may consider the possibility that only a minor part of 
the particles comes from a non-equilibrated source, theoretically this doubt may be 
extended for the whole process. Dynamical calculations indicate that the time may 
not be enough to thermalize the whole spectrum before the break up 

In this paper we would like to test a model of hadronic and quark matter, which con- 
siders a non-exponential one-particle distribution. Such distributions with power- 
law tail belong to a class of anomalous distributions, not satisfying the central limit 
theorem of statistics. In particular the canonical Tsallis distribution, (1 + E/b)~ c , 
which has been derived by an attempt to generalize the traditional thermodynam- 
ics [0], extrapolates between the expontial curve for E <C b and the power-law tail 
for E ^> b. The apparent temperature of the spectra is given by T — b/c. The 
Tsallis index is related to the power as q = 1 — 1/c or q = 1 + 1/c depending on the 
treatment of energy and particle number. 

1.3. Long tailed distributions 

Before reviewing the main formulae of the non-extensive thermodynamics and ap- 
plying it for a massless quark-gluon plasma we remind to the central limit theorem 
and show a possible exemption from it. This exemption, the Lorentz distribution, 
has a Fourier transform corresponding to exponential transverse mass spectra after 
regularization. 

Considering n independent random variables their properly scaled sum is distributed 
more and more Gaussian; in the n — > oo limit exactly Gaussian. The only require- 
ment beyond independency is that the fiducial variables must have a short range 
distribution. 

In cases of unusual distributions, the central moments may diverge, or just not 
vanish for large foldness n. This is the case for the Lorentz-distribution, w{x) = 
(T/tt)/(1+T 2 x 2 ). Its Fourier transform is given by, w(k) = exp(— \k\/T), whose log- 
arithm cannot be differentiated at k = 0. The distribution P n {x) of x — a n ^2 i=1 Xi 
has the Fourier transform P n (k) and 



log P n (k) = -n\k\a n /T. 



(1) 
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The scaling a n = 1/n leads in this case to a result independent of n and hence to a 
limiting distribution, which happens to be the same Lorentzian. 

Furthermore a physically supported regularization with a finite mass m renders the 
central moments finite. 



^)=expf^f±^j (2) 

leads to a distribution of the arithmetic sum x = (l/ n )S™=i x ii (phyiscally inter- 
pretable as a center of mass or energy), with a Fourier transform, 



P n (k) = *«(*/„) = exp ^ m -^ +(nm)2 j . (3) 

This may be the one-particle spectra of hadrons recombined from n partons each 
with an almost vanishing mass m — > 0, but with a finite total mass M = nm. A 
simple, but powerful picture of hadron formation from several partons, which leads 
to the experimentally observed exp(M — M t ) transverse mass scaling spectra. 
This cannot be, however, the whole story. Not only because of the power-law tail of 
distributions, but also on theoretical grounds. Pure recombination, as considered 
above, would need re-heating and expansion for not decreasing the entropy during 
the spontaneous hadronization. Although no detailed microdynamical description 
is yet known for this process, there are doubts, whether enough expansion may 
take place for this entropy requirement. Therefore it is of interest to circumlinc 
possible mechanisms which would change the precursor quark matter distribution 
into another one of the hadrons, closer to equilibrium and hence possesing higher 
entropy. 

Such a case has been suggested in The Tsallis distribution, 

W1 {E) = (1 + E/by c = w TS {E,T, 1 - 1/c), (4) 

after n-fold recombination becomes another Tsallis distrinution, 

w n {E) = (1 + E/nby nc = w TS (E, T, 1 - 1/nc), (5) 

much closer to the Gibbs-limit Woo(E) — wts(E,T, 1). It seems to be promising 
therefore to investigate the non-extensive thermodynamics of quark matter. 



2. Tsallis distribution and energy per particle 

In this paper we concentrate on the question that how large the energy per particle 
E/N (simply related to the entropy per particle at vanishing chemical potential) 
can become in a quark matter with Tsallis distribution. Assuming a large enough 
volume for the fireball and homogenecity, we take E/N = e/n, the corresponding 
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ration of densities. We shall compare this for the traditional relativistic Boltzmann 
gas with massless particles, for the original bag-model description of the quark- 
gluon plasma (QGP), and for a Tsallis-QGP with power-law tail distributed quasi- 
particles. We shall point out, that only the last version is able to reproduce the high 
value of E/N « 6T fitted to RHIC experiments by the massive hadronic thermal 
model. This makes the Tsallis-QGP to an interesting candidate for the pre-hadron 
quark matter. 

2.1. Non-extensive thermodynamics 

In the following we brcifly review basic formulae in the non-extensive thermody- 
namics. Tsallis has suggested an altered definition of entropy, which includes the 
Boltzmannian case as a limit. Typically it leads to a power-law tailed energy dis- 
tribution instead of the exponential Gibbs-distribution, known from equilibrium 
statistical physics. This seems to fit for relativistic heavy ion collisions, where the 
observed transverse momentum spectra of hadrons definitely show a power law tail. 

The Tsallis entropy contains a parameter q, the familiar thermodynamics is recov- 
ered in the limit q = 1. It is given by 



with Wi being the probability of the occurrence of the state i in the statistical sample 
of physical states. The parameter q can be related to a power characteristic to the 
high-energy part of these probabilities by q = 1 + 1/ 'c. The entropy can be expressed 
with the help of c as follows 



The canonical probability is obtained by maximizing the entropy J7J) with a con- 
straint on the energy and, if applies, on some charge conservation: 



Instead the sum of the energies of the individual states, E' — ^ ji WiE i , and corre- 
sponding charges N' — J+ WiQi in acertain version of Tsallis thermodynamics one 
considers the q-weighted average quantities as thermodynamical variables: 




(G) 




(7) 



S q - @{E- fiN) = max. 



(8) 



E = 



EXj 
EX 




N = 



EX 



(9) 



One assumes that these E and N are additive for subsystem and reservoir pa- 
rameterizing this way non-linearities (correlations) in energy and particle number. 
(Tsallis brings only one argument in favor of this choice against the naive one: the 
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latter "turns out to be inadequate for various reasons, including related to Levy- like 
superdiffusion, for which a diverging second moment exists" . No reference found to 
this, yet.) 

The optimal canonical distribution can so be obtained as 

Wi = ~exp c (-a:i), (10) 

with Xi = [Ei — fiQi)/T. Here we introduced the parameter T, which we refer to 
as "temperature" in the following discussion. This is, however, not equal to the 
inverse of the Lagrange multiplier /3 [ E| • 

The basic thermodynamical relation can be expressed, as follows: 

- Q = P V = T\n c Z = TS- (E-fiN)(l - S/c) (11) 

with ln c (x) = c(l — x — 1 ' c ). This relation defines the pressure, and nicely resembles 
the thermodynamics in the usual form with corrections in the order of 1/c. 

The pressure is a function of T and fi only in the homogeneous limit. The derivatives 
of the grandcanonical thermodynamical potential are 

^PV=P, ^P V = S ' ^pV = N(l-S/c), (12) 

This helps to re-express the energy and the particle number in terms of the pressure 
and its derivatives: 

* = V— ^ = V °"_v%. - ( 13 ) 

1 c dT 1 c 8T 

It is interesting to note, that the energy per particle E/N does not contain the 
C(l/c) correction factor directly, it resembles the usual expression 

E _ V-dy. + 1 dT P 

We note here, that using the E' = J] WiEi average energy for the canonical con- 
straint (and similarily N' — ^ WiQi) leads right away to the same distribution, 

vii = ~(l+Xi/c)- c , (15) 

but now q = 1 — 1/c (the relation to the Tsallis index changes). Also the thermo- 
dynamics, defining pV by T\n c Z contains finite 0(1/ c) modifications: 

TS =pV(l + S/c) + E - fiN. (16) 
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Now also the derivatives of the grandcanonical potential receive corresponding cor- 
rections, but the expressions for E and N are formally the same, as in the other 
case. The relation of the entropy to the pressure, however, changes, 

dp 

S = V 7^vM> (17) 

1 c dT 

as well as the definition of pressure now belongs to a q index reflected to 1. 

Summarizing, irrespective to which Tsallis-case we choose, we define the pressure 
and the grandcanonical potential correspondingly, so that in both cases 

P V = Tc(l- Z- 1 ' ). (18) 

The canonical distribution shows a power-law tail: 

w i = ^ r ,, 1 , v= - (19) 
Z (1 + Xi/c) 

The partition sum Z is obtained by requiring ^\ ti>j = 1. 



2.2. Tsallis ideal gas, bosons, fermions 

In the homogeneous case, very often encountered in simple physics models, the 
partition function is considered as a product of contributions from each particle 
state with a given momentum k. The familiar logarithm of such a partition function 
is then additive in the phase space. 

In the case of the non-equilibrium Tsallis-distribution the entropy is, however, no 
more additive; the above mechanism is no more trivial. Assuming that the energy, 
volume and particle number are still additive, the Tsallis ideal gas the following 
pressure: 

/rl 3 k 
J^Tln c Z, (20) 

with Z being the one-particle Tsallis partition sum, 

Z = ]T cx Pc (-nx), x = (E k - n)/T. (21) 

n 

here exp c (x) is the inverse function of ln c (x). This is the basis of the distinction 
between Fermi-Dirac and Bose-Einstein statistics: in the former case only n = 
and n = 1 contributes to the sum, while in the letter case all non-negative n-s from 
zero to infinity. 

The 1/c-uncorrected derivatives of the pressure give, 

f d 3 k gfc-M f d 3 k , Ek-ji 

n, = J e = J J^E k n c (^-) (22) 
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with 

n c (x) = -^ln c Z(x). (23) 

Let us now analyze the one-particle distribution and the partition sum. First we 
regard the Boltzmann limit, exp c (— x) <C 1. We get 

d d ( x\ ~ c ~ 1 

n c (x) w -— ln c (l + exp c (-a;)) w -— exp c (-a;) = 1 + - . (24) 
cte eta; V c/ 

It is seen, that the Tsallis temperature, T and the inverse slope of the logarithmic 
one-particle spectra somewhat differ for finite c in this scenario: 

(c+l)T slopc = cT = b. (25) 

The product quantity, b is regarded as an energy (or transverse momentum) scale 
which characterizes the onset of the power-law tail of the canonical Tsallis distri- 
bution. The general shape is given by 

n c {E) = (l + E/b)- c - 1 (26) 

for vanishing chemical potential /j = 0. In this case the particle number density 
and energy density integrals can be calculated analytically, resulting 

1 2b 3 _L 66^ 

U °~ 27r 2 c(c-l)(c-2)' 60 ~ 27r 2 c(c-l)(c-2)(c-3) { ' 

This suffices to express a relation between the inverse slope of the one-particle 
energy spectrum at small energy E <C b, the energy per particle and the Tsallis 
parameter q=l+l/c: 

E/N = e/n - T slopc . (28) 

An estimate can be given from the RHIC experiments, pointing out an inverse slope 
of 175 MeV (cleaned from transverse flow effects) and E/N = 1 GeV. This agrees 
with q = 1.135 in the Boltzmann ideal gas limit of massless partons, a much more 
likely candidate for the pre-hadron matter, then the usual QGP: at the same E/N 
value the massless Bose distribution leads to T = 325 MeV temperature and slope. 

2.3. Finite chemical potential 

We are interested in the constant E/N curve on the T — \x plane for a mixture of 
massless bosons and fermions, which are Tsallis-distributed. Denoting (E — qmu)/T 
generally by x we deal with the following type of integral experessions: 



n 



= [dkk 2 n c (x)= [ dx (xT + fi) 2 n c (x). (29) 
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Fig. 1. Energy per particle divided by the temperature (E/NT) as a function of 
the Tsallis index. 

Similarily for the quantity eo, just a higher power of (xT + /i) occurs under the 
integral sign. The one-particle distributions n c (x) are for fermions or for bosons 
respectively. 

Denoting the general integrals by I n (a) — f^° a dx x n n c (x), the result consists of 
different power contributions. For massless fermions and antifermions together we 
get 

n = T 3 S 2 (a) + 2T 2 l j,D 1 {a) + 2> 2 S (a), (30) 

and 

eo = T 4 S 3 (a) + ST 3 fiD 2 (a) + +3T 2 p 2 S 1 (a) + Tfi 3 D (a). (31) 

with S n (a) = I n (a) + I n {— a) and D n (a) = I n (a) — I n {— a), a = [i/T. For bosons 
in chemical equilibrium one takes fj, = 0, so only the leading integrals remain, of 
course, in this case with Bose-Tsallis n c (x) functions. The energy per particle for a 
mixture, like quark matter, becomes 

E_ _ d B e 0B (T) + d F e aF {T, + d F eo F (T, -fi) 
N ~ d B n 0B (T) + d F n 0F (T, /j.) + d F n 0F (T, -p)' 

For c — > oo one is back with the familiar Fermi-Dirac and Bose-Einstein distribution 
functions, and both cq and no are finite order polynomials of T and [i. The ratio 
is a growing function of both, so a constant E/N reflects in a convex curve in the 
T — /i plane. For large enough c the Tsallis version has to behave similarily, but 
this question we can only decide numerically. 

3. Comparison of Quark Matter Models 

In conclusion we have compared three models for quark matter. The simplest, the 
massless Boltzmann gas, as it is long known, gives E/N = e/n = 3T. The equation 
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Fig. 2. Energy per particle for massless Boson gas with Tsallis and Gibbs distri- 
butions as function of the spectral slope. 
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Fig. 3. The E/N ~ 1 GeV line on the temperature - chemical potential plane for 
the Gibbs and for the Tsallis distribution. 



of state is given by, 

p = \<tT 4 , e = crT 4 , p = nT. (33) 

Supporting the Stefan-Boltzmann gas with a bag constant, B, we arrive at the 
traditional bag-model QGP. Now the situation is a little bit more complex. The 
equation of state is modified, 

p = -o-T 4 - B, e = crT 4 + B, p = nT - B. (34) 
3 

This kind of matter has a stability limit at p = 0, so the bag constant is related to 
a critical temperature B = The energy per particle becomes at this limit 

E/N = e cr /n cr = AT. 
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Finally in a Tsallis-distributed quark matter the Stefan-Boltzmann constant a de- 
pends on the power c (or Tsallis index q = 1 — 1/c) . The equation of state becomes 

p=\a{c)T i , e = o-(c)T 4 , n = A(c)T 3 . (35) 

The energy per particle in the Boltzmann limit is given by 

E/N = — C —ZT 1 (36) 
c — 4 

which is about 6T for c = 8. In a bag constant supported Tsallis-QGP version one 
obtains, E/N = icT/(c - 4) and a power of c = 12. 



Acknowledgement (s) 

This work has been supported by the Hungarian National Research Fund, OTKA 
(T034269). 



References 

1. J.Letessier, J.Rafelski, A.Tounsi, Phys.Lett. B328 (1994) 499-505; G. Torrieri, 
J. Rafelski, New J. Phys. 3, 12, 2003; P. Braun-Munzinger, J. Stachel, J.P. 
Wessels, N. Xu, Phys.Lett. B344 (1995) 43-48; P. Braun-Munzinger, J. 
Stachel, J. P. Wessels, N. Xu, Phys.Lett. B365 (1996) 1-6; M. Bleicher, F. M. 
Liu, A. Kernen, J. Aichelin, S.A. Bass, F. Becattini, K. Redlich, K. Werner, 
Phys.Rev.Lett. 88 (2002) 202501; W.Broniowski, A.Baran, W.Florkowski, 
AIP Conf.Proc. 660 (2003) 185-195; 

2. J. Cleymans (Cape Town, South Africa) , H. Oeschler (TU, Darmstadt), K. 
Redlich (Wroclaw, Poland), Phys.Lett. B485 (2000) 27-31; K. Redlich, S. 
Hamieh, A. Tounsi, J.Phys.G27:413-420,2001 

3. F.M.Liu, K.Werner, J.Aichelin, M.Bleicher, H.Stoecker, J.Phys.G 30: 
s589-s594, 2004; M.Bleicher, H.Stoecker J.Phys.G. 30: s lll-sll8, 2004 

4. C.Tsallis, J. Stat. Phys. 52, 479 (1988) , |cond-mat/0312500| 
cond-mat/9903356, cond-mat/0012371; C.Tsallis, E.Brigatti, 

|cond-mat /0305606| C.Tsallis, E.P.B orges,|cond-m at/0301521 C.Tsallis, 
F.Baldovin, R.Cerbino, P.Picrobon, cond-mat/0309093, C.Anteneodo, 
C.Tsallis, cond-mat/0205314; L.G. Moyano, F.Baldovin, C.Tsallis, 

|cond-m at /03050£ip 

5. T.S.Biro, B.Muller, Phys.Lett. B578 (2004) 78-84 

6. C.Tsallis, F.Baldovin, R.Cerbino, P.Pierobon, |cond-mat /0309093 



